Mesoscopic competition of superconductivity and ferromagnetism: conductance peak 

statistics in metallic grains 
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We investigate the competition between superconductivity and ferromagnetism in chaotic ultra- 
small metallic grains in a regime where both phases can coexist. We use an effective Hamiltonian 
that combines a BCS-like pairing term and a ferromagnetic Stoner-like spin exchange term. We 
study the transport properties of the grain in the Coulomb blockade regime and identify signatures 
of the coexistence between pairing and exchange correlations in the mesoscopic fluctuations of the 
conductance peak spacings and peak heights. 



Although superconductivity and ferromagnetism com- 
pete with each other, the coexistence of both states has 
recently been reported in heavy-fermion systems [l| and 
high-T c superconductors @- This observed coexistence 
has led to enormous renewed interest in those materials. 
A main objective is to find systems in which ferromag- 
netism and superconductivity can be easily transformed 
into each other by changing an experimentally control- 
lable parameter. 

In nano-sized metallic grains a coexistence regime has 
been predicted theoretically 0, 0, S] ■ It was shown that 
the presence of spin jumps in the ground-state phase di- 
agram of a single grain is a unique signature of this co- 
existence and that an external Zeeman field can be used 
to control and tune the size of this regime [5| . 

Experimentally, the competition between supercon- 
ductivity and ferromagnetism in metallic grains has not 
yet been investigated, in part, because fabrication and 
control of nanoscopic metallic devices is very challenging. 
The first nanoscopic single-electron transistors (SETs) 
with superconducting islands were produced by breaking 
nanowires into small pieces . Discrete energy levels and 
pairing effects were observed in a single grain by mea- 
suring the tunneling conductance (for a review see 0]). 
Itinerant ferromagnetism has been studied as well in nor- 
mal metallic grains such as Cobalt Q. During the last 
decade, production and gating of metallic grains has been 
considerably improved, e.g., using chemically synthesized 
grains and electromigration [^, [To| • A very recent devel- 
opment is the observation of superconductivity in doped 
silicon, which might further simplify the production and 
control of mesoscopic superconducting devices [11 1 . 

Here we calculate the tunneling conductance in the 
Coulomb blockade regime for an ensemble of almost- 
isolated metallic grains using a rate equation ap- 
proach [HI ■ Our analysis is based on an effective Hamil- 
tonian for chaotic or disordered systems, which combines 
a superconducting BCS-like term and a ferromagnetic 
Stoner-like term originating in pairing and spin exchange 
correlations, respectively 3, 14]. By renormalizing the 
pairing coupling constant, we can derive accurate results 
using exact diagonalization of the many-body Hamilto- 



nian within a truncated band. 

We propose that signatures of the coexistence of pair- 
ing and ferromagnetic correlations can be observed in the 
mesoscopic fluctuations of the conductance peaks in su- 
perconducting grains with a sufficiently strong exchange 
interaction. We find a regime in which the peak spacing 
distribution exhibits bimodality induced by pairing cor- 
relations, while the mesoscopic fluctuations of the peak 
heights are suppressed by the presence of ferromagnetic 
correlations. Comparison of experimental data and the- 
oretical predictions for such systems might serve as an 
ideal testing ground for the interplay between the charge, 
spin and pairing channels in disordered systems. 
Model. We consider a metallic grain whose single- 
particle dynamics are chaotic and with a large dimen- 
sionless Thouless conductance <?t- In the absence of an 
orbital magnetic field, such a grain is described by an 
effective universal Hamiltonian of the form 
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The one-body term in ([I]) describes the kinetic energy 
plus confining single-particle potential. Here c\ a is the 
creation operator for an electron with spin a = ± in the 
single-particle level tk lying within a band of width gx ■ 
The second term on the r.h.s. of (JTJ) is a charging energy 
term which is constant for a grain with a fixed number 
of electrons N (C is the capacitance of the grain) . The 
third term is the pairing interaction with strength G and 
where P 1 * = ^ c l+ c l- is the pair creation operator. The 
fourth term in ([1]) is an exchange interaction expressed in 



terms of the total spin operator S = ^ s 
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(tj are Pauli matrices). The parameter J s is the exchange 
coupling constant (estimated values of J s for a variety of 
materials were tabulated in Ref. [HI]). 

The pairing interaction can only scatter time-reversed 
pairs from doubly-occupied to empty levels but does not 
affect the singly-occupied levels (referred to as "blocked" 
levels) . Thus each many-body eigenstate of the Hamilto- 
nian ([1]) \i) = \U, C; B, 7, S, M) factorizes into two parts. 
One part £} is a zero-spin eigenstate of the reduced 
BCS Hamiltonian, J2ka e kc\ cr Cka — GP^P, constructed 
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within a subset U of empty and doubly-occupied levels. 
Here Q denotes a set of quantum numbers that distinguish 
the different eigenstates within U. For example, C can be 
defined as the quantum numbers of the Slater determi- 
nant that is obtained from \U, () in the non- interacting 
limit, i.e., limc^o \U, £) = |£)- The other part of the 
eigenstate 7, S, M) is obtained by coupling the set of 
singly-occupied levels B (each of these levels has spin 
1/2) to total spin S and spin projection M. Here 7 de- 
notes a set of quantum numbers distinguishing between 
states with the same spin S (TEG]}. For a specific set 
B of b singly-occupied levels, the total spin ranges from 
S = (S = 1/2) for an even (odd) number of electrons 
to S — 6/2 with each spin value having a degeneracy of 
db(S) = ( s+b / 2 ) - (5+1+6/2) • 

The only relevant energy scale of the reduced pair- 
ing Hamiltonian (J s = 0) is A/5, where A is the bulk 
pairing gap and 5 the single-particle mean-level spacing. 
The number of levels in the grain can be truncated by 
renormalizing G such that the low-energy spectrum of the 
grain remains approximately the same. For a pickctfcncc 
spectrum, the renormalized coupling constant is given 
by G r /5 = [arcsinh ((N r + l/2)<5/A)] _1 where N r is the 
number of levels in the truncated band [lij ]. Strictly 
speaking, this holds in the absence of an exchange inter- 
action. However, since the exchange interaction affects 
only the blocked levels, the renormalization holds as long 
as the number of blocked levels is small compared with 
the total number of levels in the truncated band. 

The pairing and exchange terms in (TT]) compete. Pair- 
ing correlations tend to minimize the total spin of the 
grain while the exchange interaction favors maximal spin 
polarization. This competition leads to interesting and 
non-trivial behavior that is the focus of this work. 
Conductance. We consider a grain that is weakly cou- 
pled to leads. In the sequential tunneling regime a typ- 
ical tunneling width T satisfies T <C S,T, and we as- 
sume the charging energy to be the largest energy scale 
e 2 /2C » J s , A, (5, T. The conductance displays a series 
of sharp peaks as a function of gate voltage, each peak 
describing a tunneling event, in which the number of elec- 
trons in the dot changes from N to N + 1. In the rate 
equation approach, the conductance is determined solely 
by the many-body energies and tunneling rates T,^ in 
the left (right) lead between an eigenstate i of the N- 
electr on g rain and an eigenstate j of the (N + l)-electron 
grain 12]. These transition rates are given by 
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where we have used the notation a — {LA, £, B, 7) and in- 
troduced overall coupling strengths '. The operator 
V4-(r) = SfcV'fcMc/L creates an electron with spin pro- 
jection a at the left (right) point contact ri (r r ). There 
is only one term in the sum over single-particle states k 
that can contribute to the matrix element in (T5|). This 



tunneling level (denoted by A) is always a good quantum 
number, cither in the final state j = (a',S',M') if tun- 
neling occurs into an empty level (hence n' x = 1) or in 
the initial state i = (a, S, M) if an electron tunnels into a 
singly-occupied level (hence n\ = 1). Since the total spin 
of the grain is conserved we can apply the Wigner-Eckart 
theorem and express the matrix element in ([2]) in terms 
of a reduced matrix element which is independent of the 
magnetic quantum numbers. In this case, we find for the 
reduced transition rates 
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with r^ r = Tq T \ipx(r\ iT )\ . The first case n' x — 1 is char- 
acterized by B' = B + {A} and W =14- {A}. The state 
I^a) = R\Pn x =o\M, C) is obtained from the Slater deter- 
minant expansion of \U, by keeping only those determi- 
nants in which the level A is empty (defined by the projec- 
tor P nx =o) and removing this level (defined by the opera- 
tion R\). In the second case n\ = 1 we have B' = B— {A}, 
W = U + {A}, and \V' X ) = R\P nx =2\U',C,')- The rates in 
Eq. J3|) are used as an input for a reduced system of rate 
equations, which is obtained from Eqs. (9a, b) in Ref. [12] 
by replacing the full indices with the reduced indices 
(aS,a'S'). 

In a chaotic grain (without a time-reversal symmetry- 
breaking magnetic field), the single-particle levels and 
wave functions follow the Gaussian orthogonal ensemble 
of random matrix theory [ijj]. We study the statistics 
of the Coulomb-blockade conductance peaks at a typical 
experimental temperature of T = 0.15 Q. For each re- 
alization of the one-body Hamiltonian, we calculate the 
five lowest many-particle eigenstates of the many-body 
Hamiltonian (TTJ) using the Lanczos method for a band- 
width of N r = 8 and up to 19 electrons. Using the many- 
body energies and wave functions as input, we calculate 
the transition matrix elements ([3|) and solve the system 
of reduced rate equations. For each tunneling event, the 
maximal value of the conductance as a function of gate 
voltage is found numerically. The above procedure is re- 
peated for ~ 4000 samples to collect sufficiently good 
statistics. In the limit A = we compare with larger 
bandwidth results obtained from a closed solution of the 
rate equations [3] and find good agreement. We have 
also verified that the contribution of higher many-body 
eigenstates is negligible for T < 0. 15. 
Peak spacing statistics. The peak spacing A2 mea- 
sures the distance between two neighboring conduc- 
tance peaks. Conductance peak-spacing distributions are 
shown in Fig.[T]for different values of A and J s . In the ab- 
sence of pairing and exchange interactions (A = J s = 0), 
the distribution is bimodal because of the spin degener- 
acy of the single-particle levels. The exchange interaction 
tends to suppress this bimodality since it leads to spin 
polarization and thus mesoscopic fluctuations of the spin. 
At J s /S = 0.6, the bimodality is completely washed out. 




FIG. 1: Peak spacing distributions at T = OAS for several 
values of the pairing ratio A/5 = 0,0.25,0.5,1. Each panel 
shows results for both J 3 — (dashed histograms) and J s = 
0.65 (solid histograms). The dotted line in the top left panel 
is the analytic distribution for A = J s = and T 5. 



The most significant effect of pairing correlations on 
the peak spacing distribution is to restore bimodality. 
For an exchange coupling of J s /5 = 0.6, bimodality 
starts to reappear for A/5 — 0.5. At even stronger 
pairing A/S — 1, the peak spacing distribution is well 
separated into two parts. The left part describes the se- 
quence of even-odd-even (EOE) tunneling events, while 
the right part corresponds to odd-even-odd (OEO) tran- 
sitions. For A/S > 1 the ground state spin of an even 
(odd) grain is almost always 5 = (5 = 1/2) and meso- 
scopic fluctuations of the spin can be ignored. Increas- 
ing the exchange coupling constant merely shifts the two 
parts of the distribution closer together. 

The separation of the peak spacing distribution into 
two parts in the presence of pairing correlations can be 
understood qualitatively using the fixed-5 BCS approxi- 
mation at T = H for a picketfence spectrum. For an 
EOE sequence the first peak corresponds to the blocking 
of an additional single-particle level and the second peak 
to the removal of this blocked level by creating an ad- 
ditional Cooper pair. An expansion of the BCS ground- 
state energy in powers of 5/ A yields to leading order 
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In an OEO tunneling sequence, one first creates a Cooper 
pair and then blocks a single particle level. In this case 



A 



(OEO) 
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In Fig. [2] we show that Eqs. (J4j> and |5]) are a good 
approximation to the exact average peak spacing A 2 for 
A/S > 2. Since the pairing gap enters in Eqs. (f3]) and 
([S]) with opposite signs, the separation of the two parts 
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FIG. 2: (a) Average peak spacing A2, and (b) relative peak 
height fluctuation width o-(g max )/(/ max as a function of the 
pairing ratio A/S for J s — (open circles) and J s — 0.65 
(solid circles) at temperature of T — OAS . In panel (a) we 
compare the results with the value of A2 for a picketfence 
spectrum (solid lines) and the fixed-S BCS mean-field theory 
(dashed lines). The error bars in panel (b) are the estimated 
statistical errors using 4000 random-matrix samples. 



of the distribution becomes larger with increasing pairing 
gap. This effect is reduced in the presence of an exchange 
interaction as is also seen in Fig. [2^. The exchange inter- 
action strength J s /S is a material constant and for most 
metals is smaller or comparable to J s /5 ~ 0.6 [ljj . In 
contrast, the pairing ratio A/8 can be made larger by 
using a smaller grain. Thus the bimodality effect can 
be tuned experimentally by measuring grains of different 
sizes. We also observe in Fig. [2] that a picketfence spec- 
trum provides a good approximation for the mesoscopic 
average of A 2 . While the results presented here are for 
T = OAS, the peak spacing statistics do not change sig- 
nificantly at lower temperatures. 

Peak height statistics. Peak height distributions are 
shown in Fig. [3] [we use ln(g max /g max ) where g max is the 
average conductance peak height] . In the absence of ex- 
change and pairing interaction (A = J s — 0) the dis- 
tribution is known analytically at T <C 5 20]. The 
effect of exchange interaction on the width of the peak 
height fluctuations is somewhat similar to the effect of 
temperature. While temperature increases the number 
of levels contributing to the conductance, the exchange 
interaction brings down high-lying spin states and thus 
increases the number of states contributing to the con- 
ductance at small but finite temperatures. Consequently, 
the ratio o(g max )/g max between the standard deviation 
and average of the conductance peak height decreases in 
the presence of exchange (see Fig. [2Jd). In particular, we 
observe (see Fig. [3] for A = 0) that the exchange interac- 
tion suppresses the occurrence of small peak heights 21 1 . 
This suppression was observed in semiconductor quan- 
tum dots, where the Cooper pairing channel can be ig- 
nored and a closed solution for the conductance is avail- 
able There, the finite- temperature suppression of 
the probability of small conductance peak heights in- 
duced by exchange correlations led to a better agreement 
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FIG. 3: Conductance peak height distributions 
P(ln(g max /g max )) at T — 0.18 in the presence of pair- 
ing and exchange correlations (for notation see Fig. [TJ . We 
compare with the analytic distribution for A = J s — and 
T < 8 20] (solid lines). 

between theory and experiment. 

Pairing correlations induce a gap that pushes larger 
spin states to higher energies. Thus, we observe that 
the ratio a(g max )/g ma ^ (Fig. ^jp) and the probability of 
small conductance peak heights (Fig. [3]) increase with A 
(up to A/5 ~ 1) for J s = 0.65. For a pairing inter- 
action that is strong enough to destroy all spin polar- 
ization (e.g., A/ 5 > 1 for J s = 0.65), the peak height 
distribution becomes essentially independent of the ex- 
change interaction strength and the suppression effect of 
c r (5max)/ffmax by exchange correlations completely disap- 
pears (see Fig. [2b). 

Mesoscopic coexistence. A mesoscopic coexistence 
regime of superconductivity and ferromagnetism can be 
defined by the simultaneous occurrence of a bimodality in 
the peak spacing distribution (a signature of pairing cor- 
relations) and the suppression of peak height fluctuations 
(a signature of ferromagnetic correlations). In Figs. [1] 
and [3] such a coexistence case corresponds to A/ 5 = 0.5 
and J s /5 = 0.6. Interestingly, for a picketfence spectrum 
this case still belongs to the superconducting regime [Bj]. 
Thus, we find that mesoscopic fluctuations increase the 
parameter regime for which coexistence of pairing and 
exchange correlations can be measured. The coexistence 
regime should be directly observable in Platinum, which 
has recently been shown to be superconducting in gran- 
ular form [22| and also has a relatively lar ge e xchange 
interaction strength of J s /5 ~ 0.59 — 0.72 [15[. In ex- 
periments it might be difficult to vary the exchange in- 
teraction strength (a material constant). However, we 
expect to find similar effects at fixed exchange interac- 
tion strength by tuning an external Zccman field. 
Conclusion. In this work we calculated the transport 



properties of an ultra-small metallic grain in its Coulomb 
blockade regime and studied the mesoscopic conductance 
fluctuations in the presence of both pairing and exchange 
correlations. Of particular importance is the finding that 
experimentally accessible mesoscopic signatures of super- 
conductivity and ferromagnetism can coexist in a certain 
regime. 
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